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Abstract 

We calculate connected correlators in time dependent Gaussian or- 



c3 

r^ ' thogonal and symplectic random matrix ensembles by a diagrammatic 

' I method. We obtain averaged one-point Green's functions in the lead- 

^ ! ing order 0{N^) and wide two- level and three- level correlators in the 

Cj . first nontrivial order by summing over twisted and untwisted planer di- 
agrams. 
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Introduction Recently, Brezin and Zee pointed out the universality of wide connected correlators 
in random matrix theories |l|. They recognized its importance from the viewpoint of physics in 
disordered systems, though this universality in random matrix theories for two dimensional quantum 
gravity was already shown by Ambj0rn, Jurkiewicz and Makeenko 0. This universal feature of 
the wide connected correlator is shown by explicit calculations in quite extensive classes of random 
matrix ensembles in various ways [^ ^, §, while the well-known short distance correlator can be 
calculated in limited narrow classes [|T], Q. We have to take into account these purely mathematical 
properties of random matrix theories themselves in order to recognize the essentially universal nature 
of level statistics of physical systems. 

In this letter, we calculate the wide connected two level and three level correlators in gaussian 
orthogonal (GOE) and symplectic (GSE) ensembles with time dependence as discussed in some 
unitary ensembles |0, ||, [l^, [ll| . The diagrammatic method is known to be useful to calculate wide 
connected correlators in some unitary invariant ensembles ^ |T^ and in GOE ||^. We are going 
to calculate the explicit form of the two level and three level correlator in GOE and GSE by this 
method. The results in GOE is consistent with that already obtained by Verbaarschot, Weidenmiiller 
and Zirnbauer. We review their method for correlators in GOE and extend it to those in GSE. For 
these ensembles, we have to sum over both orientable and non-orientable diagrams. The results can 
be checked to compare those calculated by solving functional equations @, [l^ or by replica method 



H|, when time dependence is suppressed. 



Gaussian orthogonal ensemble To begin with, we explore the case of the GOE which is defined 
by the ensemble of real symmetric matrices obeying the probability distribution 



P(H) = -e.p\^--J_JtTr[[^-j +m^/f1|. (1) 

where H is an N x N matrix for GOE. The measure DH is explicitly written as 

N 

DH=l[dHkkY[dH.j. (2) 

k=l i<j 

A dressed line is defined by 

which is independent of time t because of the transrational invariance of the distribution function 
P{H). The averaged one-point Green's function G{z) is calculated from this one-point function 

Gij[z) 

1 1 1 ^ 

G{^)^J^(^^T^) = J^Y.Gu{z). (4) 

4 = 1 



The one-point function Gij{z) is expanded in the power series of H 



1 OO 1 

G^A^) = -j:-MHnj (5) 



Z n=0 ^ 

We evaluate this function with the diagrammatic decomposition by using the free propagator of the 
real symmetric matrix 

{Hij{s)Hki{t)) = ^{SuSjk + Sik5ji)K{s - t), (6) 

where K{t) = 2^^"™'*'" ^^ ^^^ extend the time dependence of the probability distribution as the 
time dependent factor K{t) = cr^e""^*-* with an arbitrary suitable function u{t). The corresponding 
diagram to this free propagator which has a twisted part is shown in fig 1. The one-point function 
can be calculated by summing over all oriented planer diagrams in eq(^ . Any diagram with twisted 
propagators does not contribute to the leading order in the 1/N expansion, and therefore the same 
diagrams contribute to the leading order as in the unitary ensemble. When the series eq@ is written 
in 

G.W = ^Eh^ 3n, (7) 



^ n=0 \ / 

Qn represents the number of the oriented planer diagrams with n free propagators. The number gn 
obeys the following recursion relation 



n-l 



9n+l = 9n+J2 9mgn-m, Tl > 1, (8) 

where go = 1- The corresponding diagrammatic representation to the recursion relation @ is 
depicted by fig 2. On the other hand, the vacuum polarization function 

/ 2 \ « 
E{z)^z-G{zr'=zY.^Sni^—j . (9) 

satisfies the equation 

g{z){z-j:{z)) = i. (10) 

This gives the relation among the expansion coefficients gn and s„ 

n 
Y^ Sngn-m = 0, n > 1, (11) 

m=0 

where the coefficient Sq = —1. We can show the following relation from the eq(|ll|) by the mathe- 
matical inductivity 

n-l 
■Sn ^ fl'n+1 / J 9m9n—in- V-'-^J 



m=0 



This equation indicates the vacuum polarization function consists of the one-particle irreducible 
diagrams. This equation (|T^) and the recursion relation (H) give 

Sn = gn, n> 1, 

and therefore 

S(^) = ^G{z) = ^' (13) 

The vacuum polarization function T,{z) is obtained by solving the quadratic equation 

1 



S(z) = (z-Vz^-2a^). 

Then we obtain the averaged one-point Green's function 

Giz) = ^iz-Vz^-2a^). (14) 

The obtained function G{z) satisfies the boundary condition G{z) — > 1/z as 2; — > oo. 
Now we turn to the connected two level correlator 

G(zi, Z2; t) = ( — Tr ^—— — Tr ^— r ) • (15) 

We calculate this function in the following expression of the resolvent operator where we have to 



only sum over simpler diagrams following Brezin and Zee |12 



"1 00 00 00 "1 1 

G{z,, z,- t) = ^^7^ E E ^(Tr/7(t)-Tri7(0)"),. (16) 

N-^ ozi 0Z2 ^i ~^ mn zY'zq; 

First, we ignore contractions within the same trace, in which we take into account only m = n parts 

00 1 

n=l ^ l^l^2j 

In the nth order, n untwisted planer diagrams and n twisted planer diagrams in fig 3 contribute to 
the leading order 

1 d' ^2n(K{t)y_ 2 d' ( Kjt) \ 

N^ dzidz2 t'l ^' V2^i-S2y N^ dz,dz2 °^ 1^ ^^i^J ' ^ ' 

Next we include contractions within the same trace in (Trif (t)™Trilf(0)"). The following operator 
product expansion (OPE) formula 

(Trlog(l-if(t)M)i^.,(0)) = 2:^((^-— ^y ), (19) 

tells us that the two level correlator is obtained by replacing the free one point Green's function Sij/z 
to the dressed one Gij{z) in eq(|l^) 

2 92 / K(t) \ 

G{z,, Z2; t) = -^^^^ log (^1 - ^Giz^)Giz2)j (20) 

This result with respect to the connected part agrees with that obtained by solving functional equa- 



tions [§, and by a replica method [|14 



Gaussian symplectic ensemble The Gaussian symplectic enseinble(GSE) is the ensemble of 

Ta(i 
ij 



quaternion real Hermitian matrices. The component if" of 2N x 2N quotation real Hermitian 



matrix H is written as 

3 
ra/3 _ o-(O) ra/3 , . V- TT{a)al3 



Hf = <5"^ + z E <o-r , (21) 

a=l 

where a^ {k = 1, 2, 3; a, P = 1, 2) is a component of the 2x2 Pauli matrix, H^, is a real symmetric 
matrix and Hj^j (i, j = 1, ■ ■ ■ , N) is real antisymmetric matrix in H^j . The trace of H^ in eq.p^) 
reads 

N N 3 

TtH' = E HfHfr = 2EEE(^i;V- (22) 

«ja/3 i=l i=l a=0 

Here we define the distribution function of GSE 



^(^'4->^(-t/->"''|(^)^+"^^^i 



(23) 



with the measure 

DH=l[dH^^'l[I[dHS^- (24) 

i<j fe=l i<j 

The free propagator of H is 

{Hff{s)H2f{t)) = j^K{s - t) {6uS,,6-^6^^ + 6,,6,,A-^'''^) , (25) 

where the tensor A"^'"^^ is defined by 

^a/3,7. ^ ^a/3^7e _ §'^^§^1 _ (26) 

The first and second terms in the free propagator are represented by an untwisted and twisted 
diagrams shown in fig 4 (a) and (b), respectively The calculation of the dressed line is done by the 
expansion 

The one-point Green's function is defined as its trace 

As in the GOE case, only untwisted planer diagrams contribute the one-point Green's function in 
the large N leading order 



Since the number Qn of the oriented planer diagrams with n free propagators is identical to that in 
the GOE case, the vacuum polarization function Il(z) is calculated by summing over one-particle 
irreducible diagrams. 



E{z) = z-G{z 



,-1 



The vacuum polarization function Ti{z) is obtained by solving the quadratic equation 

S(z) = hz-Vz^-8a^). 

Then we obtain the averaged one-point Green's function 

G{z) = ^{z- v/i^^8^) . (31) 

Note that G{z) tends to 1/z as z ^ oo. 

Next we compute the two-point function G{zi, Z2) defined by 

/ 1 11 1 \ 

G{zi, Z2] t) = ( 77^7 Tr — — - 77^7 Tr- 



2N zi-H(t)2N Z2-H{0)/^ 

EE ;^{TTH{tr TTH{Oy%. (32) 



4A^2 Q-,^ Qz^ ^^ ^^ ^^ ^m.^n 

Ignoring contractions within the same trace in the following expression 



E (33) 



ra=l 



It is the same as in the GOE case that in the nth order n untwisted planer diagrams and n twisted 
planer diagrams in fig 3(a) (b) contribute to the leading order. A part of the twisted diagram shown 
in fig 3(c) is calculated by the following contraction formula for the tensor A°'^'"''^ defined by eq (p6D 



Z ^nNctictn+ij/SiAi+i — y^ Aaia2,l3il32 Aa20i3,/32l33 _ _ _ Aa„a„+i,l3„/3„+i 

a2,---,a„,l32,---,l3„ 

= 2"""'" 74°-'-°"+-'-'^-'-'^"+-'- ('34') 

This gives a trace formula 

TrA" = ^(A")""'/3/3 _ 2^^ (35) 

a,/3 

which enables us to obtain the two level correlator in this truncation 

' '' £^f^)"p".T...).-^^,oJ.-?£«). (3a) 



4A^2 dzidz2 ~| n"^ \Z1Z2 1 2N^ dzidz2 \ -21-22 



iFiom this expression, we obtain the exact two level correlator by replacing the free line 5^5°"^ / z to 
the dressed one G°'j{z) in eq(|36D to incorporate contractions within the same trace 

G{zi, z,- t) = ~^^^^hg{l-2K{t)G{z,)G{z,)). (37) 

This result with respect to the connected part agrees with that obtained by solving functional equa- 



tions [|, |1^ and by replica method |I4 



Three level correlator We can also calculate the connected three level correlator both in GOE 
and GSE. That in the GOE case is defined by 

G{zi, Z2, 2:3; ii, i2, h) = ( ^Tr — — — Tr — — — Tr — — )^ 

N zi- H{ti) N Z2- H{t2) N 2:3 - if (ts) 

1 93 ~ - ~ 1 ^^ fH{h)Y' fH{t2)Y' fH{ts)Y'\ , ^ 



N'^ dzidz2dz'i ^^^^^^nxn^U'i \ Zx ) \ ^2 / V ^3 , 

For GSE, we have to replace A^ — > 2A^ in the right hand side of the definition (pHf). First we discuss 
the GOE case. Ignoring the contraction within the same trace, a diagrams with mi, m2 and m-^ 
propagators is depicted in fig 5 and fig 6. The number n^ of operators in eq(^) can be expressed in 
terms of the number vtia of propagators in a diagram in fig 5 or fig 6 

?7,i=m2-|-m3, 712 = ma + mi, tt-s = mi + m2. (39) 

The summation can be calculated in terms of the number mi, m2 and vrtj, of the propagators. In the 
leading order of the 1/A^ expansion, there are 4^1712^3 diagrams with rrii, m2 and 7713 propagators 
which consist of nin2n^ untwisted diagrams in fig 5 and 3?7.in2n3 twisted diagrams in fig 6. Therefore 
the three level correlator is calculated as 

- - - 1 ^^ fH{t,)Y' (H{t2)Y" (H{h)Y\ 

EEE;— -(Tr^ ^M^ ^M . ^^ 

„i=l „2=1 „3=1 ^1^2'T'3 V ^1 / V ^2 / V ^3 / 



£ - - ~ (K{t2-h)Y' (K{t,-h)Y" (Kih-t^'^"^' 



^2 

^ ^I'T^i «it^i ™t^i V 22:2^3 / V 22:3^1 / V 22:12:2 



mi=l m2=l mz=\- 

4 / - - (K{t^-t^)Y' lK{t^-t2)Y' , \ <.. 

+ 7? E E ^r ^r + 2 (^Vf^^'^f^ permutaUons . (40) 

^ \m2=l m3=l V 22:32:1 ) \ 2ZiZ2 ) j 

To obtain the exact 3 level correlator, we have to take into account the vertex correction as well as 
the self energy correction to the one point function. The vertex corrections can be done with the 
untwisted propagators as depicted in fig 7. Then the three level correlator becomes 

4 d'^ 

G{zi, Z2, Z3; ti, t2, ts) = ~T^ o — o ^ — F{Zi, Z2, Z3; ti, t2, h) , (41) 

6 



where 

-^^23 -'^Sl -^12 



F{zi, Z2, Z3; ti, t2, ts) = Tj 7^-j T^Z 7^ 

i — A23 i — A31 i — A12 

X31 X12 1 X12 X23 1 X23 X31 1 

1 — -^^31 1 — -^12 1 — -'^ll 1 — ^12 1 — -^23 1 ~ -^22 1 — -^23 1 — ^31 1 ~ -^33 

with Xab = K{ta—th) I {2zaZh) in this approximation. If we take into account the self-energy correction, 
the free hne 5ij/ Za should be replaced by the dressed line Gij{za). Therefore we have to use 

Xab = ^Kita - h)G{Za)G{Zb), (43) 



where G{z) is given in eq(p!4D. The exact three level correlator in GOE is given by eq(^6D with the 
function (|42|) . This result agrees with that obtained by Verbaarschot, Weidenmiiller and Zirnbauer 



Next we calculate the three level correlator in GSE. All the contributing diagrams in this case are 
identical to those in the GOE case. The contribution from the nin2^3 untwisted diagrams becomes 

£ ~ ~ ~ 1 / 2K{t2-t,) Y' / 2K{t,-h) Y' / 2K{h-h) Y' 

nJ^,^,^^,2[ Z2Z, ) [ Z,Z, ) [ Z,Z2 I 

+ 1^ H Ho ^ + 2 cyclic permutations . (44) 

^\m2=lm3=l2\ ^3^1 / V Z1Z2 ) j 

Note that the number of loops with respect to the spin index in an untwisted diagram with 1711,1712 
and m^ propagators shown in fig 5 (a = 1, 2) is always mi + m2 + m^ — 1. Contribution from 
remaining 3ni?T,2^3 twisted diagrams in fig 6 is 

3 ^ ^ ^ 1 f 2K{t2-t,) Y' ( K{t,-ti) Y^ , ^^^,,,^^, ( K{t,-t2) Y' , ,^3,,,,,, 
+ l[g g /^ ^fe-ti) \ " ^^„,,^^a,p^f K{ti-t2) \ \j^msyp,ee ^ 2 cycHc permutattons] 

X \m2=lm:i=l\ ^^Zl J \ Z1Z2 J J 

+ ^ I E E f ^^^'~^'^ r' (A-2)-,/5/3 f 2K{ti-t2) Y' ^ 2 c^c/^c permutations] . (45) 

^\m2=lm3=l\ ^32:1 / V ^1^2 / J 

The contraction formula ( |34D for the tensor A°'^''^'^ enables us to calculate these contributions, and 
then the three level correlator in GSE is obtained as 

1 d^ 

G{zi, Z2, Z3; ti, t2, ts) = -TTTI o a a -^(^1' ^2, Z3; ti, ts, h) ■ (46) 

AN^ OZ1OZ2OZ3 
The function F{zi, Z2, z^; ti, t2, t^) is defined by eq(^) with the function Xab for GSE 

Xab = 2K{ta-tb)G{Za)G{Zb), 



with the one point Green function defined by eq(|3lD. Three level correlators in GOE and GSE 
obtained by the diagrammatic method are identical those calculated by the replica method and the 



functional method 13 



We have calculated the averaged one point Green's functions and the wide connected two level 
and three level correlators in Gaussian orthogonal and symplectic random matrix ensembles with 
time evolution by a diagrammatic method. All results are consistent with those obtained in all other 



methods M, |, 13, 14 
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Figure 1: The diagram of free propagator in the GOE case, (a) an untwisted free propagator, (b) a 
twisted one. 
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Figure 3: (a) An untwisted planer diagram in the nth order with the definition of n untwisted 
contraction parts, (b) A twisted one with definition of n twisted contraction parts. In the case of 
the GOE, the spinor indices {ai,Pi,. . . ) should be neglected. 
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Figure 4: The diagram of free propagator in the GSE case, (a) an untwisted free propagator, (b) a 
twisted one. 





Figure 5: An untwisted diagram with mi,m2 and 7713 propagators, (a) mi,m2 and ^37^0 case, (b) 
mi = 0, 1712 and ^7237^0 case. 
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Figure 6: A twisted diagram with 1711,7712 and 7713 propagators, (a) 7711,7112 and ^7237^0 case. (b),(c) 
and (d) mi = 0, 7712 and ^7137^0 case. 
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Figure 7: Diagrams which require vertex correction. (n>l) 
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